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ABSTRACT 

We study steady state winds from compact objects in the regime where the wind velocity at infinity is 
ultra-relativistic. This may have relevance to some models of Gamma-Ray-Bursts (GRB's). Particular 
attention is paid to the case where neutrinos provide the heating. Unless the neutrino luminosity is 
very large, L > 10 54 erg/s, the only allowed steady state solutions are those where energy deposition is 
dominated by neutrino-antineutrino annihilation at the sonic point. In this case, the matter temperature 

£S) ' near the neutron star surface is low, less than lMeV for typical neutrino luminosities. This is in contrast 

to the case for sub-relativistic winds discussed in the context of supernovae where the matter temperature 
near the neutron star approximates the temperature characterizing the neutrinos. We also investigate 
the setting of the neutron to proton ratio (n/p) in these winds and find that only for large (> lOMeV) 
electron neutrino or electron anti-neutrino temperatures is n/p entirely determined by neutrino capture. 

OO ' Otherwise, n/p retains an imprint of conditions in the neutron star. 

Subject headings: Stars: Neutron - Gamma-rays: Bursts - Neutrinos 

(N ■ 

J> ' 1. INTRODUCTION 

ly-j ' We examine properties of ultrarelativistic winds from compact objects and we thereby obtain insights into the special 
, case where these winds arise from neutrino heating. Our principal motivation for the present study is to determine 
^C) ■ if physically plausible steady state solutions exist in the ultrarelativistic regime and to determine the relation between 
C^) ' neutrino parameters and the character of the outflow. A peculiar feature of the steady state ultrarelativistic solutions is 
t-H , that the matter temperature near the neutron star surface is cold. We comment below on whether or not this is physically 
■ realizable. We also study how the final neutron to proton ratio (or alternatively electron fraction, Y e ~ l/(n/p+ 1)), is 
(— i | set in these models. Only for large (> lOMeV) electron neutrino or electron anti-neutrino temperatures does Y e in the 
O !~ initial wind material come into steady state equilibrium with the neutrino fluxes. 

There is a rich literature concerning steady state outflow from neutron stars. In connection with supernovae the 
neutrino-driven wind occuring several seconds post core bounce is interesting because it is the favored candidate for the 
site of the r-process elements (Meyer et al. 1992; Woosley et. al. 1994; Takahashi, Witti, & Janka 1994). The relation 
between neutrino parameters and the wind have been extensively discussed (Duncan, Shapiro, & Wassermann 1986; Qian 
& Woosley 1996). These studies were concerned with gentle, modest entropy winds, characterized by Lorentz factors 
of order unity far from the star. Near the neutron star surface these gentle winds are well approximated as being in 
hydrostatic equilibrium. The temperature is set by the competition between neutrino heating and neutrino cooling terms, 
i.e., the matter temperature is roughly equal to the temperature characterizing the neutrinos. Apart from effects due to 
cd nucleosynthesis, the final electron fraction in these gentle winds is set by the ratio of electron neutrino (u e ) and electron 
anti-neutrino (v e ) capture rates on nucleons (Qian et al. 1993). The general relativistic extension of the wind equations, 
redshift factors, and neutrino trajectory bending have been discussed by Fuller & Qian (1996), Cardall & Fuller (1997), 
and Salmonson & Wilson (1999). These works also emphasize subrelativistic outflow velocities. Detailed studies of 
relativistic effects along these lines for r-process nucleosynthesis have been carried out recently by Otsuki et al. (2000). 

In a pioneering study Paczyhski (1990) examined ultrarelativistic How from neutron stars relevant for studies of GRB's. 
In this work no reference was made to a heating term and, except for photon diffusion which is unimportant near the 
neutron star surface, the flow was assumed adiabatic. This work provided a detailed picture of the wind evolution far 
from the neutron star surface, which is the regime relevant for optical emission in GRB's. There are many similarities 
between these steady state ultra-relativistic winds and the fireballs discussed in connection with GRB's (Shemi & Piran 
1990). 

The present study is an attempt to fill the gap between detailed studies of the relation between neutrino heating 
and sub-relativistic outflow and the study of ultra-relativistic outflow without reference to a heating mechanism. This 
is also a detailed look at the 'guts' of a simple model for a GRB central engine. It may be applicable to the model 
proposed by Salmonson, Wilson, & Mathews (1999). In this model binary neutron stars are compressed and heated as 
gravitational radiation saps orbital kinetic energy during the last seconds of inspiral. The heating drives copious neutrino 
loss (L„ ~ 10 52 ergs _1 ) which in turn drives a hot pair plasma from the neutron star surface. 
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Lastly, it has recently been argued that the electron fraction in neutrino-heated fireballs may be a clue to neutrino 
physics at the central engine, and may have observable consequences (Fuller, Pruet, & Abazajian 2000). The present 
study serves as a "proof of concept" of this idea for a simple case. Interestingly, we find that the electron fraction in 
the wind far from the central engine is set by processes very near the central engine. If this is a generic feature of GRB 
central engine models, information about the electron fraction in GRB's may help pin down the central engine. 

2. EQUATIONS 

The general relativistic equations describing the outflow of material from the neutron star are Euler's equations: 



u • (VT) = -u • (VT„) 
(g + u ® u) • (VT) = - (g + u ® u) ■ (VT„) , 



baryon number conservation, 



V • (p b u) = 0, 

and the change in electron fraction due to lepton capture on baryons, 



u 



(VY e ) — (1 — Y e ) (A„ e n— > p + A e + n _> p ) — Y e (Ap eP _> n + A e -p_> n ) . 



(1) 



(2) 



(3) 



->p < "e+n— >pj *e \'*v e p— »n i '»e _ p- 

In the above u is the 4- velocity of the outflowing material, T and T„ are the stress-energy tensors for the outflowing 
matter and the neutrinos respectively, g is the metric tensor, and p\, is the rest mass energy density of baryons in the plasma 
comoving frame. In Eq. (3) A„ en -> P and A PcP ^ n are the capture rates for electron anti-neutrinos and electron neutrinos on 
neutrons and protons respectively, and A e + n ^ p and A e - p ^ n are the capture rates for positrons and electrons on neutrons 
and protons respectively. In what follows we employ the Schwarzschild metric and use natural units, G = c = h = kb = l, 
with fcfe the Boltzmann constant. In Schwarzschild coordinates the line element is 

ds 2 = g^dx^dx" = - (1 - 2M/r) dt 2 + (1 - 2M/r)~ 1 dr 2 + r 2 dtf, (4) 



Here M is the mass of the neutron star and dfl 2 = d6 2 + sin 2 



Imposing steady state conditions and writing the 

.1/2 



radial component of the 4-velocity u r as u r = vy where y = (1 — 2M/r) 1 ' 2 / (l — v 2 y" and v is the outflow velocity of 
the plasma as measured in an inertial frame at rest in the Schwarzschild coordinate (t, r, 9, <j>) system allows us to rewrite 
the equations above as 



(1 - v 2 )ap h = (p + P) 



M 



1-v 2 r 2 (1 - 2M/r) 



QPb = vy 



P' + (P + P) 



M 



+ P' 



1 



v(l - v 2 ) r 2 (1 - 2M/r) 



M = 4irr p^vy 



(5a) 

(5b) 
(5c) 



vyY' e = (1 - Y e ) (A„ eIwp + A e+Iwp ) - Y e (A PeP ^ n + A e - p ^ n ) . (5d) 

In the above, primes denote differentiation with respect to r, and p and P are the total energy density (including 
rest mass) and pressure, respectively, of the plasma in the comoving plasma frame. We have also introduced a = 
— ((g + u <g> u) ■ (VT„)) r /pb and q = — u • (VT„) /pb which are respectively the momentum and energy deposition into 
the plasma per unit time per unit baryon rest mass. M is the mass outflow rate. The second law of thermodynamics and 
some algebra allows us to produce from these expressions two more first integrals of the flow: 

ds m b q 



dr Tvy 
d{Hy) 



(6) 



; =y((l-v 2 )a + q/(vy)). (7) 

Here we have introduced the nucleon mass Wb, the entropy per baryon s (in units of Boltzmann's constant), and the 
enthalpy H = (p + P)/pb of the plasma. It is useful to have explicit formulas for these quantities: 



Ts 

H-1 = H = — (8b) 

In this equation p$ is the density in units of 10 8 gcm~ 3 . The number 5.21 above becomes 1.89 as T decreases below the 
threshold for pair production. In writing Eqs. 8a and 8b we have assumed that the electrons are non-degenerate and that 
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the entropy of light particles dominates the entropy of nucleons (i.e. s > 20). This is relevant for our analysis because we 
will be concerned principally with conditions near the sonic point, and we will show below that low entropy conditions 
near the sonic point cannot lead to ultrarelativistic flow at infinity. The quantity H is related to the quantity 77, the ratio 
of energy density to rest mass energy density often discussed in connection with GRB's, through 77 ~ (3/4) ff + 1. 

If the flow is to be relativistic at infinity it must either be supersonic at the neutron star surface or pass through a 
sonic point somewhere above the neutron star. In the present work we focus on the case where the flow is not already 
supersonic at the neutron star surface. For the case of nondegenerate electrons and negligible baryon pressure (s > 20) 
the equation governing v' can be written as 



{ (l + H)(v 2 -v 2 s )v'= l - 



%H-(\ + 2H/3) (M/r) + ar (l - v 2 ) - 

3 v 1 — 2M/r v ; 3u r 



v (1 - v 2 ) K sj r 

From this equation we see that at the sonic point the following conditions must be satisfied, 



(9) 



y2 = v s = YW+T ^ 



(1 - .9 (r)) \H = g(r) + ^--(l-v 2 ) ar. (10b) 



and 

C\-a(r\) . . , 

3 av ' 3u 
In this last equation we have introduced 

g (r) = Mj (r - 2M) (11) 

and v s 2 = {dp /dp) is the sound speed with the differentiation is taken at constant entropy. 

Eq. (10b) implies that at the sonic point qr/u r H < g(r). This is interesting because the change in entropy of the 
outflowing gas past the sonic point is implicitly given by 

f ds f qm h dr f qdr ( qr \ 

J T = J T7^ = J U^ff < Iconic ~ ^ (12) 

Here the subscript s denotes values at the sonic point and the second to last inequality holds if q/u r H decreases more 
rapidly than r _1 . When the inequality in Eq. (12) holds, then, the entropy is approximately constant past the sonic 
point. Because the hydrodynamic equations imply that 3T'/T = s'/s — 2/r — (u r '/u r ), the inequalities in Eq. (12) also 
imply that the temperature cannot increase substantially past the sonic point for accelerating flow. This implies that 
the enthalpy per baryon is roughly constant or decreasing past the sonic point. We then arrive at the important result 
that that the final Lorentz factor of the wind is not much larger than the value of H at the sonic point when the above 
inequality is satisfied. (The final Lorentz factor of the wind occurs when when the energy in relativistic particles ~ Hp^ 
has been converted to baryon kinetic energy.) This will allow us to make important conclusions regarding the sonic point 
conditions needed in order that the final flow be relativistic. 

We can roughly divide the possible sonic point solutions into two families: those for which the terms involving the heating 
and momentum deposition terms q and a may be neglected (Family I), and those for which these terms are important 
in determining the conditions at the sonic point (Family II). As we will show that the neutrino energy deposition terms 
satisfy the inequalities in Eq. (12) when H and u r increase with radius, only family II solutions can support steady state 
ultrarelativistic flow. An exception to this occurs when the sonic radius is very near 3M. Previous discussions of steady 
state winds from neutron stars have focused on family I solutions. For these solutions H ~ (3/2)M/ (r — 3M) at the sonic 
point, and hence the sonic point must be near 3M in order to produce ultrarelativistic flow in the absence of heating 
terms (Paczyhski 1990). For neutrino driven wind studies relevant for the r-process the heating terms are generally not 
important near the sonic point simply because those flows are subrelativistic or mildly relativistic, so that the sonic point 
occurs far from the neutron star where neutrino heating has fallen off. 

For Family II solutions (those where heating determines the sonic point conditions and drives the flow relativistic) 
some general observations can be made. If we assume that H(r s ) 3> H(ro) and s(r son j c ) ^> s(ro) (where r son i c is the 
schwarzschild radial coordinate of the sonic point and r is the initial radial coordinate of the flow, usually taken to be 
the neutron star radius here) then the sonic point condition and Eq. 7 imply 

r sonic yq _ r Bonic q ( yqr \ n 

In Eq. 13 we have defined y e ff, we have assumed that (1 — g)H ^> 1 at the sonic point, and we have neglected the 
momentum deposition term a. Eq. 13 serves to specify the position of the sonic radius given a heating term q and a 
velocity profile. If q/u r oc r _/3 (which occurs in steady state winds if the heating rate per unit volume drops as r -/5 ~ 2 ) 
the sonic point occurs at 

2(0- Iconic , (14) 



r \2/e//(l - ffsonic) 
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which implies that the sonic point occurs near Tq for rapidly dropping heating terms. Again the momentum deposition 
term a has been neglected in this approximation. Another interesting property of these winds is that the temperature 
does not vary greatly within the sonic radius if q is a decreasing function of radius. This is seen by noting that Eqs. 6 
and 13 imply 

\ f r "" ic g/Tu r 



which allows us to to define an effective temperature within the sonic point, T e ff = T s y s /y e ff. 

The above equations allow us to extract the properties of ultrarelativistic steady state winds for arbitrary heating 
functions. We now turn to the particular case where neutrinos from the central compact object are responsible for the 
heating. 

3. NEUTRINO DRIVEN ULTRA- RELATIVISTIC WINDS 

Perhaps the first question raised by a study of steady state neutrino-heated wind solutions is whether or not these 
winds are physically realizable for conventional compact objects. There are two parts to this question: i) whether or not 
steady state can be achieved; and ii) whether or not the flow can be relativistic. We defer consideration of the second part 
until later, when we discuss mass ablation from the neutron star surface. While a hydrodynamic simulation is needed 
to determine the timescale to achieve steady state in these winds, a rough estimate is that the equilibration timescale 
is of order the sound crossing time between the sonic radius and the neutron star surface, ~ 10~ 4 s. This timescale is 
short compared to the neutrino diffusion timescale, which is of the order of seconds. The neutrino diffusion timescale 
governs the evolutionary timescale of the neutron star unless it undergoes a phase transition or becomes dynamically 
unstable. It is known that for most cold neutron star equations of state general relativistic instability sets in for r ~ 3M. 
Below this radius the star becomes dynamically unstable and collapses on a dynamic timescale which is comparable to 
the equilibration timescale given above. For these small radius stars, then, it may not be sensible to put much stock in a 
steady state wind solution. 

Before describing in some detail properties of the allowed steady state solutions, we must first examine the neutrino 
energy deposition terms. To parametrize the neutrinos we make the usual assumption that the neutrinos are emitted 
from a neutrinosphcrc with some radius tq (in this work we do not make a distinction between neutrinosphcre radius and 
neutron star radius) and are characterized at the neutrinosphcre by a Fcrmi-Dirac black body energy distribution with 
temperature T and zero chemical potential, 

This is a crude approximation to the actual expected energy spectra but suffices to mock-up the energy deposition 
physics (see Cardall & Fuller (1997)). We further make the free streaming approximation, i.e. we neglect the fact that 
the neutrinos suffer a small depletion with increasing radius owing to interactions with the outflowing plasma. The facts 
that real neutron stars are characterized by some small but finite decoupling region over which the neutrino distribution 
function continues to evolve and that neutrinos from real neutron stars are in general non-thermal or better characterized 
by degenerate spectra are not important for our arguments. For a discussion of when it is apropriate to treat the neutrino 
as sharply decoupling from a neutrinosphere see Salmonson & Wilson (1999). 

Evaluating VT„ for various neutrino plasma interactions is most easily accomplished by noting that Liouville's the- 
orem implies that in an inertial frame at rest at some Schwarzschild coordinate r the neutrino distribution function is 
characterized by the red-shifted temperature 

T v (r) = T u (r ) ( \~ 2 ^° ) ' = hT v {r ), (17) 



1 - 2M/r ) 

(see Fuller & Qian (1996)) which also serves to define the redshift factor h, and a maximum angle of deviation from the 
radial direction given by 

T \ 2 1 - 2M/r \ 1/2 



r 



cos max = x= 1 - I — — — . (18) 



1 - 2M/r J 



The mass M of the neutron star appears in the expresion for x because of the bending of null trajectories in the 
Schwarzschild geometry as has been discussed by Cardall & Fuller (1997) and Salmonson & Wilson (1999). 

In an inertial frame at rest in Schwarzschild coordinates, the contribution to VT„ = (q3p,Q£p) from ncutrino- 

antineutrino annihilation is 

Q * = 9^ CTv (r)4 T " (r)4 Fa (0) Fs (0) ^ (r) + Tp (r)) $ {X) 

Q »> = W^f CT " M 4 T * M 4 F, (0) F 3 (0) (T„ (r) + T v (r)) * (*) /4 (19) 
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In the above $(ar) = (x - l) 4 (x 2 + 4x + 5), = (x - l) A {x + l)(3x 2 + 9x + 8), G F is the Fermi constant, and 

hats are used to denote quantities in an inertial frame at rest in the Schwarzschild coordinate sytem. The quantities 
i*3(0) = 77r 4 /120 and ^(O) « 23.3 are the Fermi-integrals of argument zero. The quantity C depends on the neutrino 
species. For v e v e annihilation C = (1 + 2 sin 2 8 w ) 2 + 4 sin 4 6 W where 6 W is the Weinberg angle (sin 2 9 W « 0.23), while 
for VpVp and v T v T annihilation C = (1 — 2 sin 2 W ) 2 + 4 sin 4 9 W . For simplicity in what follows we will assume that the 
luminosity of a single neutrino flavor gives all vv annihilation heating and that for this flavor the neutrino and antineutrino 
temperatures are equal, T v = T P . We will refer to this single neutrino flavor as fj, or r. This is merely a calculational device 
designed to reproduce the total heating rate from all neutrino flavors and adjusted accordingly for sensitivity studies. 

Making a transformation back to Schwarzschild coordinates, we find that pb<3W = (1 — v 2 )~ 1 / 2 Q vj? (l — w^E'/4$) and 
PbO-vv = (1 — u 2 ) _1 (l — 2M/r)~ 1/ ' 2 (2°p(5 , /4<f> — v). The terms entering the right hand side of the sonic point condition 
are then 



r«„e _ KT 9 r 6 (T u (r)/MeV) 9 ^ H ( 



3u r (1 - v 2 y/ 2 u r (T/MeV) 4 V 4 * 

, 2 , 3-10- 8 r 6 (T„(r) /MeY)\ ^ { \ 

^ (1 V ] T = (l-2M/r)V 2 (T/MeV ) 4 * H U 'I ' (20) 

where r 6 = r/10 6 cm. 

Evaluation of the contribution to (VT„) from neutrino-plasma interactions is complicated by the fact that an iner- 
tial observer comoving with the outflowing plasma sees a neutrino distribution function whose temperature varies with 
direction as T„(r)(l — vcos9)/(\ — v 2 ) 1 / 2 . The effect of this velocity (and direction) dependent rcdshift is to make the 
neutrino-plasma heating rates small as the outflow velocity becomes large. We will show that neutrino-antineutrino anni- 
hilation is the dominant heating term for the conditions of interest here. It suffices for us to note that in the limit where 
v = 0, the contributions to (VT„) from neutrino-electron scattering and neurino capture on baryons are, respectively, 

T 4 I 1 — x 

^~ 10 <7-7^-^- (21) 

if T < T v , and 

T 4 (TjY e +Tl(l-Y e )) l-x 
**~ 10 *" 5- (22) 

Note that all of these heating terms satisfy the inequality in Eq. (12) if H increases with radius. Neutrino- Antineutrino 
annihilation satisfies Eq. (12) because q ui >/H ~ $r 2 /(M_ff) which falls more rapidly than r~ 6 . Neutrino electron 
scattering satisfies Eq. (12) because q ve j 'H <~ (1 — x) <~ r~ 2 Neutrino capture on baryons also satisfies this condition 
because q v b/H ~ (1 — x)/H. These estimates for the scaling of the heating rates with radius neglect the dependence of 
the heating on the redshift factor h. This implies that the heating rates decrease even more rapidly with increasing radius 
if the wind is accelerating. This means that if the final flow is to be ultra-relativistic, we must have H > 1 at the sonic 
point, i.e., either r son i c ~ 3M, or, 

T ~ 1 MeV ^^d-^y/lOM.V)^ ^ ^ 
for vv annihilation dominating the heating, or 

/I - \ ^ 5 

T v > 40 MeV h- 1 - 9 - J (23b) 



r 6 {l - x) 

if electron-neutrino scattering dominates the heating rate, or, 

T v > 30 MeV h- 1 I — 9 — ) (23c) 

\r 6 {l-x)J 

if neutrino capture on baryons dominates. In Eqs. (23a, 23b, and 23c) the vvalues of all quantities are taken at the sonic 
point. 

Unless the sonic radius occurs near 3M (g ~ 1), the neutrino temperatures needed to drive steady state ultrarelativistic 
flow imply extraordinary neutrino luminosities. (For example, consider the cases inherent in Eqs. 23b and 23c, L w 
10 54 r 2 (T„/20MeV) 4 erg/s.) In what follows we therefore restrict ourselves to the case where vv annihilation dominates 
the heating at the sonic point and where only smaller, more realistic, neutrino luminosities are required. 



6 



3.1. Steady state flow when vv annihilation dominates at the sonic point 

Having determined that for modest neutrino luminosities (L < 10 54 ergss _1 ) ultrarelativistic steady state flow only 
occurs if vv annihilation heating dominates at the sonic point, we now turn to a discussion of the behavior of these 
solutions. 

For vv annihilation the heating scales as q/u r oc $r 2 /i 9 + 0(v) and the sonic radius is readily estimated. This is shown 
in Fig. 1. From this figure we see that the sonic radius is generically near ro- Because q„v/u r decreases rapidly with 
radius we may approximate y e ff ~ y(fo) and T e ff T(ro) « (y SO nic/yef f)T son ic, which implies that the temperature 
is "cold" within the sonic radius so long as thermal neutrino energy losses are unimportant. Here "cold" means smaller 
than the neutrino energy distribution temperatures. 

The mass ablation rate for vv dominated solutions is 

M « 10- (M 0S -) (^(l-. 9 )(^)/10McVn (24) 

V / sonic 

If M is too large the flow cannot be relativistic. This is the "baryon loading problem" discussed in connection with 
GRB's. The relation between M and the neutrino heating rates depends in detail on the structure of the compact object 
atmosphere below, above, and through the neutrinosphere. As argued by Fryer & Woosley (1998), obtaining a sufficiently 
low M may be particularly difficult for vv energy deposition because this heating rate drops so rapidly with radius. This 
requires the density scale height above the neutrinosphere to be small. For type II supernovae, calculations show that 
M « E v /(GM/ro), where E v is the integrated neutrino energy deposition rate (Woosley et. al. 1994). This implies that 
most of the neutrino energy deposition goes into extracting baryons from the gravitational potential well of the neutron 
star. For supernovae, then, ultra-relativistic winds are not expected unless the early-time density scale height is smaller 
than expected or the late time neutrino luminosities are larger than expected. Fryer & Woosley (1998) have argued that 
the baryon loading problem cannot be overcome for mass ejection from strange stars resulting from the phase transition of 
a cooling or accreting neutron star while Salmonson, Wilson, & Mathews (1999) have argued that neutrino heating during 
the compression of an inspiralling neutron star might lead to ultrarelativistic flow. The question must be considered for 
each proposed GRB site. 

An estimate of T e ff allows a simple determination of the influence of neutrino capture on the electron fraction in the 
outflow. Noting that the number of neutrino captures per baryon n c is approximately given by dn c — dr q^ / {T v \x r ) , 
where v = v e or v e for neutrino capture on neutrons or protons, gives the total number of neutrino captures per baryon 
between the neutron star surface and the sonic point 

n c « l(r 6 r 6 (fc 6 (l - g)(l - x)) sonic (Tl + T«) In ( Ssonic / S o) , (25) 

where s son j c and s are the entropy at the sonic point and neutron star surface, respectively. With this expression we 
can solve for the final electron fraction in the fireball, neglecting e + /e _ capture because of the low plasma temperature. 
Representive solutions are shown in Fig. 2. This figure serves to illustrate the basic features of how the final electron 
fraction is set: i) For low T„ c , T„ e , Y e simply remains what it was in the neutron star, ii) As T Ve , T Pe increase, the final Y e 
tends toward equilibrium with respect to neutrino capture (i.e. Y e — > T% /(T^ +T^J), while still retaining some memory 
of the value of the electron fraction in the neutron star crust. And iii) For T Vt , or T Pe > lOMcV, Y e is set ultimately by 
the competition between v e and v e captures. If the neutrino temperatures are high and neutrinos set the electron fraction 
in the outflow, the final electron fraction is expected to be low because the neutrinos emitted from a neutron rich star 
generically satisfy the temperature hierarchy T Ve < T„ e . Because neutron star crusts typically have Y e < 0.1, the electron 
fraction in the outflow is also expected to be low even when neutrino capture does not set the final electron fraction. 

A numerical example of a steady state neutrino heated wind leading to ultrarelativistic flow is shown in Fig 3. For this 
example we have taken M = 1.4M Q , r = AM, T V<1 = T P(i = lOMeV, and n = A/3H = 100 at the sonic point. As the 
wind behavior shown in figure 1 is quite independent of H (except of course for the baryon density which scales trivially 
with H) as long as H > a few, we have only shown a single numerical example. The sonic radius can be chosen to match 
a particular entropy at the neutron star surface. Increasing r son i c increases the entropy at the surface. The value of the 
surface entropy is very sensitive to changes in r SO nic- Numerically, for r son ic too large, the scale height for velocity changes 
becomes very small and the integration fails. Physically, the sonic radius is set by the requirement that the energy flow 
E at the sonic radius equals the net neutrino energy deposition rate interior to the sonic radius. Our plot only extends 
out to r w 3r . Past this radius neutrino energy deposition is unimportant and the flow satisfies the simple scaling laws 
given in Piran, Shemi, & Narayan (1993). The final Lorentz factor of the wind is approximately 150. 

4. CONCLUSIONS 

We have extended the study of steady state winds from compact objects to include the case where the wind is driven by 
a heating term and the wind velocity at infinity is ultra-relativistic. We find that for heating rates per unit volume which 
drop more rapidly than r~ 3 , the sonic point condition is dominated by the heating and momentum deposition terms, the 
sonic radius is near the compact object radius, and the temperature within the sonic radius is roughly constant. 

Particular attention has been paid to the case where the wind is driven principally by neutrino energy deposition. 
For these winds neutrino anti-neutrino annihilation dominates the heating unless the neutrino luminosity is very large 
(~ 10 54 ergs _1 ). Interestingly, then, for a given neutrino luminosity there are an infinite family of steady state supersonic 
solutions. On one branch of this family is the usual subrelativistic solution discussed for supernovae, and on the other 
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their exists a continuum of ultra-relativistic solutions dominated by neutrino heating at the sonic point. Of course, the low 
mass ablation rate needed for the ultra-relativistic solutions may not be compatible with the structure of the atmosphere 
of the compact object near the neutrinosphere. This needs to be investigated for each compact object. 

The way in which the electron fraction is set in the steady state winds we have discussed depends sensitively on 
the neutrino luminosity, with the number of neutrino captures per baryon varying as the fifth power of the neutrino 
temperature and depending only logarithmically on the final Lorentz factor of the outflow. For high v e or v e temperatures 
the final electron fraction is set by neutrino capture. Otherwise Y e remembers its value in the neutron star crust. At least 
for a simple case, then, the electron fraction in the outflow is a diagnostic of conditions within the central engine. This 
is interesting because a priori one might guess that the electron fraction comes to equilibrium at Y e = 1/2, or is always 
dominated by neutrino capture, as is the case for supernovae. It is also interesting because recently (Bahcall & Meszaros 
2000) it has been argued that neutrinos arising from pion decay in GRB fireballs may be detectable. This neutrino signal 
depends, along with the other parameters characterizing the fireball, on the electron fraction in the fireball. If more 
realistic GRB central engine models also leave their finger prints in the electron fraction in the outflow, such a signal 
might be used to distinguish between central engine models. 
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Fig. 1. — Approximate position of the the sonic point as a function of ro/M for neutrino anti- neutrino annihilation dominating 
the heating term within the sonic radius. The momentum depostion term a has been neglected and it is assumed that at the sonic 
point H(l - g) » 1. 




Fig. 2. — Final electron fraction in the fireball as a function of T„ e and for different ratios Tp e /T„ e (labeled next to the curves). 
The initial electron fraction is taken to be 0.01, and the entropy of the plasma is assumed to increase by a factor of 10 as the 
plasma travels from the neutron star surface to the sonic point. 
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Fig. 3. — Evolution of a neutrino heated wind solution as described in the text. The solid curve is for temperature (in MeV), 
the short dashed curve for q vi >, the long dashed curve for U r , and the dot-dashed curve for E = MHy. The asterices correspond 
to the position of the sonic point. 



